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Pr ■ Abstract 

Tj- . We perform a lattice study of heavy baryons, containing one (Af,) or two 

^ ! 6-quarks (5;,). Using the quenched approximation we obtain for the mass of 

Afe 

Ma, = 5.728 ± 0.144 ± 0.018GeV. 

O . The mass splitting between the A^ and the B-meson is found to increase 

by about 20% if the light quark mass is varied from the chiral limit to the 
strange quark mass. 
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1 Introduction 

Considerable progress has been achieved in the computation of low lying 
hadronic states from lattice QCD (LQCD). It turned out that quenched 
LQCD reproduces the hadron spectrum in the light quark sector surpris- 
ingly well, once finite volume corrections are carefully taken into account [||]. 
This holds in particular for the ratio mj^ /nip, which has been a notorious 
problem for LQCD over quite some years. 

Beauty physics is attracting much attention because the origin of CP 
violation in the B-system is still an open question. The investigation of such 
a system is a great challenge to LQCD. Considerable progress has been made 
in the lattice studies of heavy-light mesons like D- and B- mesons . 

Heavy-light systems from the baryonic sector so far have been studied 
with lattice methods exploratively in the limit of infinite heavy quark mass 

[|, i- 

Throughout the current study the mass of the heavy quark has been kept 

finite. We will present results for the mass of A^, a baryon composed of a 6- 
quark and two light quarks, as well as for the mass of S^, a baryon composed 
of two 6-quarks and one light quark. 

As in the case of heavy-light mesons, the most dangerous source of sys- 
tematic error stems from the fact that one is forced — in order to enter the 
region of heavy quark masses — to actually push the heavy mass close to the 
very limit of current lattice resolutions. We will attempt, however, to keep 
control on these effects of a finite lattice spacing a by a variety of precautions: 

1. We avoid to compute the masses of the A{, and the H^ directly, but 
rather calculate the mass splittings Aa = M\^ — Mb and A^ = M^^^ — 2Mb, 
with respect to the B-meson mass Mb- These splittings do not depend on 
the heavy quark mass in the infinite mass limit and are therefore less prone 
to contamination by finite a effects in the b and c quark mass regions. 

2. We monitor the dependence of the splittings on the lattice spacing 
for our three (3 values, /? = 5.74, 6.0, 6.26. This enables us (a) to check the 
assumption of weak a dependence and (b) to perform an a — i> extrapolation. 

3. We will not calculate directly the mass splittings too close to the 
b quark mass. Instead we stop the calculation at approximately twice the 
charm quark mass and then extrapolate our data to the b mass. 

Moreover, we investigate finite size effects at /5 = 5.74, on three lattices 
of spatial volumes 8^, 10^ and 12^, in lattice units. Our lattice parameters 



are detailed in table 1. To obtain a good signal for the ground state, we 
use smeared gauge invariant interpolating quark fields [Q, defined for the 
standard Wilson action. 



Ns 


Nt 


no. configs. 


/? 


arup 


8 

10 

12 


24 
24 
24 


175 
213 
113 


5.74 
5.74 
5.74 


0.542 ±0.014 


12 


36 


204 


6.00 


0.355 ±0.016 


18 


48 


67 


6.26 


0.260 ±0.014 



Table 1: Lattice parameters: space and time extents Ns and Nt, number of 
configurations, f3 and the p-mass in lattice units. 



2 Smearing Techniques and Volume Effects 

For the reasons given above, we will compute the baryonic masses with refer- 
ence to the mass of the B-meson, by proper combinations that would elimi- 
nate the 6-quark mass in the heavy quark limit. So we consider the splittings 
Aa = Ma^ —Mb and A= = M^^ —2Mb, respectively, in the single and double 
beauty sectors. In the first case, we need to compute the correlators, for the 
heavy-light pseudoscalar meson. 



Cp{t) = Y.((Qi^hV{x)) (q'ioh'Qio) 



(1) 



[X, 



and for the A baryon Q] 

X 

where Q{x) = Q{x,t) is the heavy quark field, and q"^ 



.X 



(2) 
q'^{x,t) is the 



light one, to which smearing of type J has been applied^. C is the charge 

^The correlator for the S is obtained from the one for the A, by replacing the light 
quark gjf by a heavy one Qt- 



conjugation operator given by C = i'-j'^'-)'^. 

Given the lattice results for these correlators, we perform a direct fit to 
their ratio 

RK{t) = §^ - A^-"-'' (3) 

in the large t limit. 

It is by now well known that smearing 0, 0, § is crucial to obtain a decent 
overlap of the operators with the ground state. In this work, we make use of 
the experience acquired previously while computing properties of the heavy- 
light pseudoscalar states [§], where we found gauge-invariant 'Gaussian type' 
wave functions (of r.m.s radius 0.3 fm) to provide sufficient overlap. The 
smearing was applied to the light quark source in the mesonic case. In this 
study, we are using precisely this procedure for the heavy-light baryons as 
well, without any further optimization. 

In order to establish ground state dominance we look for a plateau in 
the local mass of the ratio -Ra. In fig. la we show as an example the local 
mass of i?A, for the two largest lattices: the solid line shows the fitted value 
for the plateau. Fig. lb shows the corresponding local mass for the ratio 
R'Eif) = C^it) / Cp(t) used to determine A=. These figures show the quality 
of the plateaus for representative k values. Worse quality is found only in 
a few cases, and it resulted in larger statistical errors; for instance the two 
largest k values at /? = 5.74 given in table 3. 

The pseudoscalar mass was extensively studied in ref. and the values 
have been taken from there. 

For checking the finite volume effects, we computed Aa on three lattices 
with Nt = 24, /? = 5.74 and sizes Ng =8, 10, 12, for m = 0.156, and Kh = 
0.125, 0.140, 0.150. We compare the results for the splitting in table |^. The 
values exhibit deviations of at most 4% . In fig. 2 we plot the dimensionless 
ratio A/^/rrip as a function of L, the lattice size in units of nip] as can be seen, 
for these values of k^, the finite size effects of this ratio are smaller than our 
statistical errors. 

In the following, we will fix the volume to about 1 fm (which corresponds 
to Ns —8, 12 and 18 for (3 = 5.74, 6.0 and 6.26 respectively) and carry out a 
detailed study of the extrapolation to the continuum limit. As a tribute to 
possible finite size effects we will add the above maximal variation of 4% as 
an uncertainty to all results. 



m 


f^h 


Aa, 8^ X 24 


Aa, 10^ X 24 


Aa, 12^ X 24 


0.156 


0.125 


0.564 ±0.008 


0.554 ±0.005 


0.569 ±0.005 




0.140 


0.592 ±0.007 


0.575 ±0.005 


0.588 ±0.004 




0.150 


0.621 ±0.007 


0.601 ±0.005 


0.612 ±0.006 



Table 2: Results for three different spatial lattice extents, Ns 

Nt = 2A, (3 = 5.74. 



=8, 10, and 12, 



3 Continuum Limit 



A Splitting. The results for Aa at the three f5 values 5.74, 6.00 and 6.26 at 
fixed volume are listed in table 3, all in lattice units. The extrapolations to 
u and s-type light quarks, — for given heavy quark k/^ — are performed as 
described in 0. 

Since we are evaluating masses, it is natural to use rup to set the lattice 
scale. The values of rup for the various lattices have been listed in table 1; 
the experimental value used is: rrip = 768 MeV. 

In fig. 3a we plot Aa in GeV extrapolated to the chiral limit as a function 
of 1/Mp in GeV~^. Within the statistical precision achieved in this compu- 
tation, the points show no dependence on a, although a is varied by about a 
factor two (cf. Table 1). 

In the continuum, the 1/Mp expansion for Aa gives 



Ar'iMp)=co + ^ + 0{Mp 



-2^ 



(4) 



Assuming no dependence on the lattice spacing, this form can be fitted 
directly to the points in fig. 3a yielding Aa^ = 431(28) MeV at the mass of 
the B-meson. This result is included in the figure as the inverted triangle. 
The error bar of this point does not account for the fact that the simulation 
results exclude an a-dependence only within their precision. 

A realistic error that includes the uncertainty of extrapolating the lattice 
data to the continuum is obtained by allowing for the leading a-dependence 
at each value of Mp [^]: We start out from a selected value of Mp (in physical 
units) and interpolate the lattice results from each (fixed) /3- value to the value 
of Mp. This enables us to compare Aa at different values of a. A subsequent 



X 24, P = 5.74 



t^h 


oAa 


aAAirUs) 


aAs 


aAsirUs) 


AA(m,)/AA(m„) 


0.06 


0.25 ±0.13 


0.37 ±0.08 


-0.24 ±0.02 


-0.31 ±0.02 


1.46 ±0.40 


0.09 


0.34 ±0.10 


0.42 ±0.06 


-0.17 ±0.02 


-0.23 ±0.02 


1.23 ±0.18 


0.125 


0.38 ±0.04 


0.46 ±0.02 


-0.12 ±0.02 


-0.17 ±0.01 


1.20 ±0.07 


0.14 


0.42 ±0.02 


0.50 ±0.01 


-0.01 ±0.01 


-0.009 ±0.009 


1.18 ±0.04 


0.15 


0.47 ±0.04 


0.54 ±0.02 


0.09 ±0.01 


-0.080 ±0.008 


1.13 ±0.05 



12^ X 36, (3 = 6.00 



K^h 


oAa 


aA/^{ms) 


aAs 


aAs{ms) 


AA^rUs) / AAirUu) 


0.1 


0.22 ±0.02 


0.28 ±0.01 


-0.14 ±0.02 


-0.18 ±0.01 


1.23 ±0.05 


0.115 


0.21 ±0.02 


0.28 ±0.01 


-0.10 ±0.01 


-0.147 ±0.008 


1.29 ±0.08 


0.125 


0.25 ±0.01 


0.300 ±0.008 


-0.09 ±0.01 


-0.130 ±0.007 


1.19 ±0.03 


0.135 


0.27 ±0.01 


0.310 ±0.007 


-0.06 ±0.01 


-0.101 ±0.006 


1.17±0.03 


0.145 


0.31 ±0.02 


0.348 ±0.007 


0.01 ±0.01 


-0.049 ±0.009 


1.11 ±0.04 







18^ X 48, (3 = 6.26 




t^h 


oAa 


aAA(ms) 


cAh 


aAsirris) 


AA(m,)/AA(m„) 


0.09 


0.15 ±0.01 


0.188 ±0.009 


-0.07 ±0.01 


-0.107 ±0.008 


1.25 ±0.07 


0.10 


0.15 ±0.02 


0.19 ±0.01 


-0.073 ±0.008 


-0.110 ±0.006 


1.29 ±0.07 


0.12 


0.16 ±0.01 


0.201 ±0.007 


-0.062 ±0.06 


-0.010 ±0.004 


1.25 ±0.04 


0.135 


0.18 ±0.01 


0.216 ±0.008 


-0.032 ±0.005 


-0.071 ±0.004 


1.21 ±0.04 


0.145 


0.21 ±0.01 


0.243 ± 0.009 


0.022 ±0.006 


-0.024 ±0.005 


1.14 ±0.04 



Table 3: The mass splittings Aa and Ah at the chiral hmit and at the strange 
quark mass in lattice units for f3 = 5.74, 6.00, 6.26. The last column gives 
the ratio AA{Tns) / Aa{itLu) for the same f3 values. 



linear extrapolation in a will then yield the continuum estimate for Aa, at 
the chosen physical value of Mp. 

Reiteration of the procedure on a set of masses Mp determines the nu- 
merical dependence of A™"^*' on Mp. In order to remain in compliance with 
the assumption of linearity in a, we used only a subset of the data listed in 
table 3, by excluding the data for the two heaviest masses at any value of (3. 

A final fit of the continuum limit values A™"*' to eq.@ yields 

Aa^ = 458 ± 144 ± 18 MeV . (5) 

Here the first error is purely statistical and stems from a full jacknife analysis 
of our data. The second error represents the 4% uncertainty due to finite vol- 
ume effects. Together with the known experimental value Mb = 5.27 GeV, 
we obtain an estimate for the mass of A;, : 

Ma, = 5.728 ±0.144 ±0.018 GeV. (6) 

We note in passing, that the parameter cq is in agreement with the estimates 
obtained directly in the static approximation 0, albeit within the larger 
uncertainties of the latter. 

In the same way we obtain the mass splitting Aa,(^«s) where the light 
quark mass is extrapolated to the strange quark mass. The data fitted are 
shown in fig. 3b, and the result of the extrapolations described above is 

Aa, (m,) = 597 ± 91 ± 24 MeV , (7) 

with the same meaning of errors. Instead of looking directly at the splitting 
at the strange quark mass, one may consider the ratio 

^ AA(m,) 

AA(mO' ^""^ 

which is expected to have an even weaker a-dependence. Indeed the ratio for 
the different j3 values shown in fig. 3c follows a "universal curve" with quite 
small statistical errors. A linear fit to the (weak) mass dependence yields 

r = 1.25 ±0.03 ±0.05. (9) 

This ratio represents the change in the A?, mass when replacing the u quark 
by an s quark. 

7 



S Splitting. The mass splitting for tlie H baryon is investigated using the 
same techniques as for the A. Table 3 displays the results after extrapolation 
of the light quark to the d and s mass respectively. We then convert all data to 
physical units, and plot Ac vs. 1/Mp at the chiral limit in fig. 4. In this case 
the data show a statistically significant a-dependence, especially for heavy 
masses. We extrapolate to the continuum limit as above for pseudoscalar 
masses Mp below the charm mass. We use the same ansatz as in the case of 
Aa to fit the 1/Mp dependence of Ah and obtain: 

Ah, = -90 ± 170 MeV , (10) 

This mass shift together with the experimental value for Mb determines 

Mh, = 10.45 ±0.17 GeV , (11) 

for the physical Sfc-mass. 

Another possibility to compute the S^ mass, is to consider the ratio 
R 's{t) = C-={t)/C^{t), which should yield the splitting between the S;, 
mass and the T-meson mass. We evaluated this quantity but the resulting 
plateaus turned out to be of worse quality than those for the other chan- 
nel studied; this is due to the fact that smearing was applied only to light 
quarks, and so the heavy-heavy channels suffer from relevant contamination 
by higher-mass states. 

4 Discussion 

The mass splitting technique is a viable method to compute the A;, mass on 
the lattice. Both the lattice spacing dependence of the mass splitting and 
its dependence on the heavy quark mass are weak. Thus an extrapolation 
to the continuum and to the b-quark mass is possible. Our actual value 
of Ma^ = 5728 ± 144 ± 18 MeV can be compared to the value 5630 MeV 
suggested by Martin et al. within the naive quark model approach [Q. 

Experimentally, the determination of the A^ mass has a somewhat con- 
troversial history ever since 1981 |]12|. Recently, the UAl collaboration has 



measured the mass from 16 events in the decay channel J/\E'A|. Their value 

^It is an open question, why this decay channel has not been observed in the CDF 
experiment. 



IS 



m Ma, = 5640 ± 50 ± 30MeV. The DELPHI collaboration is presently 



quoting [|1J] a preliminary mass value, Ma^ = 5635_29 ± 4MeV, which is 



based on one candidate event in the A+tt^ and in the D'^pir— decay modes. 
These numbers are in agreement with our result. 

For M\^ we obtain 2433 ±88 ±18 MeV, which is in rough accord with the 
experimental value of 2285 MeV |]12[. Looking at the light quark dependence 



we found a 20% increase as we lift the quark mass from the chiral limit to the 
strange quark mass. Finally, an estimate for the Sf, mass is given in eq.(ll). 
Our errors do include - as the dominant part - the uncertainty induced 
by an extrapolation to the continuum limit. Nevertheless, it is desirable to 
further check these extrapolations through simulations with higher lattice 
resolutions and/or different lattice actions. 
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Figure captions 

Figure la 

The Aa local masses, given by fi^^^{t) = log[ ^ (.^ ] for ki = 0.1525, 
Kh = 0.125 at /9 = 6.0, for a 12^ x 36 lattice, and for ki = 0.1492, 
Kh = 0.12 at /5 = 6.26, for a 18^ x 48 lattice. 

Figure lb 

The same as figure la, but for /i='^(t) = log[ |,„(^n ]. 



Figure 2 

A/i^/rUp is shown vs L for three heavy quark masses aX (3 = 5.74. The 
light quark mass was fixed to about twice the strange quark mass. 

Figure 3a 

The A mass splitting is shown vs 1/Mp in GeV^^ at the chiral limit. 
The solid line is a global fit assuming no a effects. It gives the value 
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shown with the inverted triangle at the B-meson mass. The value ob- 
tained after extrapolation to the continuum limit and to the B-meson 
mass, is shown by the open circle. This extrapolation is discussed in 
the text. 

Figure 3b 

As for figure 3a but at the strange quark mass. 

Figure 3c 

The ratio A\{ms) / ^hiTny) is shown vs 1/Mp. The notation is the 
same as for figure 3a. 

Figure 4 

The S mass splitting, A= = M-^ — 2Mp is shown vs 1/Mp in GeV^"*^ at 
the chiral limit. The notation is the same as for figure 3a. 
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